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Abstract 


A  density  profile  vs  altitude  from  95  to  130  km,  obtained  from  drag  measure¬ 
ments  made  directly  from  within  a  falling  sphere,  are  presented.  The  main  sensors 
are  linear  accelerometers  aligned  to  measure  X,  Y,  Z  components  of  acceleration 
of  the  sphere  from  which  total  drag  is  obtained.  The  method  consists  of  ejecting  an 
inflatable  sphere  somewhere  between  80  and  100  km  on  the  ascent  of  am  Aerobee 
150  rocket  flight.  From  ejection,  the  sphere  continues  to  approximately  250  km 
altitude  on  a  trajectory  which  is,  except  for  drag,  a  free-fall  trajectory. 

The  prediction  of  sphere  drag  coefficients  with  their  estimated  probable  errors 
are  given  in  the  altitude  range  from  90  to  160  km  and  was  made  by  Prof.  Frederick 
S.  Sherman  of  the  University  of  California.  The  telemetry  record  showing  the  drag 
forces  used  in  calculating  the  density  profile  in  the  altitude  range  95  to  130  km  is 
shown  in  Figure  10.  The  sources  of  error  in  the  calculated  density  are  analyzed  and 
tabulated.  The  calculated  density  in  kilograms  per  cubic  meters  vs  altitude  in  kilo¬ 
meters  is  compared  with  the  ARDC  Model  Atmospheres  of  1956,  1959,  and  1961. 
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Upper-Atmosphere  Density  Obtained  From 
Falling  Sphere  Drag  Measurements 


1.  INTRODUCTION 

This  paper  gives  the  density  profile  vs  altitude  from  95  to  130  km  obtained  from 
drag  measurements  made  directly  from  within  a  falling  sphere.  The  method  consists 
of  ejecting  an  inflatable  sphere,  2.74  m  in  diameter,  somewhere  between  80  and 
100  km  during  the  ascent  of  an  Aerobee  150  rocket  flight.  At  the  center  of  the  sphere, 
supported  by  an  inflated  cylincrical  strut  (see  Figure  1),  are  the  measuring  instru¬ 
ments  and  the  necessary  electronics  to  telemeter  the  information  continuously 
throughout  the  flight.  The  main  sensors  are  linear  accelerometers  aligned  to  measure 
the  X,  Y,  Z  components  of  acceleration  of  the  sphere  from  which  total  drag  is  ob¬ 
tained.  From  ejection  altitude,  the  sphere  continues  to  an  altitude  of  approximately 
250  km  on  a  trajectory  which  is,  except  for  drag,  a  free-fall  trajectory.  Measure¬ 
ments  obtained  by  this  method  are  supplemented  by  tracking  systems  for  velocity 
and  position. 

The  drag  acceleration  aD  is  related  to  the  c  ag  force  D  and,  in  turn,  the  drag 
force  is  a  function  of  the  air  density  by  the  classical  equation 
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Figure  1 .  Inflated  Sphere  with  Instrument  Package 
Suspended  at  Center 

in  which  m  is  the  mass  accelerated  by  the  drag  force,  A  is  the  sphere  cross- 
sectional  area,  V  is  the  sphere  velocity,  is  the  coefficient  of  drag,  and  p  is  the 
atmospheric  density.  In  this  experiment,  m  and  A  are  known,  a^  and  V  are  meas¬ 
ured,  and  CD  or  p  can  therefore  be  calculated  if  one  is  assumed  or  determined 
independently.  The  authors  believe  at  this  time  that  the  drag  coefficients  available 
are  more  accurate  than  the  present  estimates  of  atmospheric  density  in  the  altitude 
range  investigated. 

This  rocket  firing  is  the  second  of  a  series  designed,  constructed,  and  instru¬ 
mented  by  project  personnel  of  the  Geophysics  Research  Directorate,  Air  Force 
Cambridge  Research  Laboratories. 

For  economical  reasons  the  first  flight  was  a  prototype  and  was  successfully 
flown  in  a  Nike  Cajun  sounding  rocket  (AA6.  170)  on  23  February  1961.  The  flight 
information  and  experimental  results  of  the  first  flight  of  the  falling  sphere  system 
for  measuring  upper-air  density  was  reported  by  Faucher,  Procunier,  and  Stark.  * 

Previous  uses  of  the  falling  sphere  were  developed  by  the  University  of  Michigan. 
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The  first  approach  was  to  eject  an  inflatable  sphere,  4  ft  in  diameter,  from  an 
Aerobee  rocket  and  measure  the  trajectory  by  means  of  the  Dovap  tracking  system4 
and  a  Dovap  transponder  in  the  sphere.  This  method  is  limited  by  the  accuracy  of 
tracking,  which  is  within  a  few  feet. 
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The  second  approach  was  to  measure  acceleration  directly  from  within  the 
sphere,  using  a  transit-time  accelerometer  in  an  aluminum  sphere  measuring  7 
inches  in  diameter  and  weighing  11  pounds.  The  accelerometer  compares  the  posi¬ 
tion  of  the  drag-accelerated  sphere  with  the  position  of  a  free-falling  mass  within 
the  sphere.  This  is  in  contrast  to  measuring  the  total  acceleration  by  a  ground 
tracking  device  and  subtracting  the  acceleration  of  gravity  to  find  the  small  difference, 
namely,  the  drag  acceleration.  This  second  approach  results  in  improved  sensitivity, 
even  though  a  small  (7 -inch)  metal  sphere  is  used  in  place  of  the  4-ft  inflatable 

o  e  e 

sphere  which  has  a  more  favorable  mass-to-area  ratio.  '  ’ 

To  utilize  this  technique  at  significantly  higher  altitudes,  it  was  necessary  to 
increase  accelerometer  sensitivity  and/or  decrease  the  mass-to-area  ratio  of  the 
sphere  as  well  as  to  propel  it  to  a  higher  velocity.  Single-axis  accelerometers 
were  used  because  they  can  provide  more  information.  The  transit-time  accelero¬ 
meter,  which  gave  average  acceleration  over  an  elapsed  time,  may  give  the  correct 
result  but  only  if  spurious  accelerations  are  negligible.  As  may  be  seen  in  Figure  10, 
page  IS ,  this  is  not  necessary  in  the  case  of  a  continuous -reading  instrument. 


2.  JNSmUMKKTATlOH 

The  test  vehicle  for  this  experiment  was  an  Aerobee  150  sounding  rocket, 

AB3.  344.  It  carried  a  payload  consisting  of  a  folded  instrumented  sphere,  radar 

beacons,  telemetry  transmitter,  a  de-spin,  and  range  safety  extensions.  The  sphere 
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was  constructed  of  1 -mil  Mylar  laminated  with  2. 25  oz/ yd  dacron  scrim.  When 
released,  it  expanded  to  a  diameter  of  2. 74  m.  The  inflated  sphere  had  an  internal 
pressure  of  3  mmHg  with  100  cc  of  methanol,  and  a  pressure  of  100  mmHg  was 
obtained  in  each  half  of  the  strut  with  24  cc  of  Isopentane.  Inflation  time  of  the 
sphere  and  strut  was  approximately  8  sec.  Figure  1  shows  the  inflated  sphere  with 
its  instrument  package  suspended  at  center. 

The  folded  sphere  was  ejected  from  the  rocket  by  a  'slingshot1  arrangement, 
utilizing  four  strips  of  neoprene  that  were  stretched  across  the  diameter  of  the  nose 
cone  (see  Figure  2).  The  sphere  was  ejected  with  a  force  of  approximately  300  lb, 
giving  an  ejection  velocity  of  approximately  9  ft/sec. 

Figure  3  shows  the  instrument  package.  It  consists  of  an  aluminum  cylinder  with 
caps  and  collapsible  turnstile  antennae,  and  is  6  inches  in  diameter,  20  inches  long, 
and  IS  lb  in  weight.  Figure  4  is  a  block  diagram  of  the  instrumentation,  which  in¬ 
cludes  three  separate  subcarrier  channels  of  orthogonal  acceleration  measurements 
and  three  time-shared  accelerometers  plus  two  internal  temperature  measurements 
on  a  fourth  subcarrier  frequency.  The  main  sensors  are  six  Donne r  Scientific 
Company  Model  4310  linear  accelerometers.  They  operate  as  subminiature  servo 


Figure  2.  Inflatable  Sphere  Container  and  Nose 
Cone 


Figure  3.  Sphere's  Instrument  Package 
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systems  responsive  to  Input  linear  acceleration  along  their  sensitive  axes. 

Basically  the  accelerometer  is  comprised  of  a  seismic  system  and  associated 
position-error  detector,  a  restoring  mechanism,  and  a  servo  error-signal  amplifier. 

The  rocket  and  sphere  were  tracked  by  use  of  two  AN/FPS-16  and  one  AN/MPS-19 
radars;  and  one  AN/FPS-16  was  also  u-ed  in  skin  tracking  the  sphere  after  ejection. 

In  addition  nine  contravies  phototheodolites  were  used  to  track  the  vehicle  and  the 
sphere,  while  two  MATTS  (Multiple  Airborne  Target  Trajectory  System)  stations 
tracked  the  telemetry  signals  from  the  sphere. 

The  accelerometers  are  calibrated  by  tipping  the  units  into  the  earth's  gravity 
field.  The  local  vertical  determines  the  maximum  value  of  one  'G'  or  9. 8035  m/sec^, 
depending  upon  local  value.  As  the  sensitive  axis  is  lowered  from  this  peak,  de¬ 
creasing  values  of  acceleration  appear  as  inputs  to  the  unit  in  the  relationship 
A  *  G  sin  0,  where  A  is  the  acceleration  input  to  the  accelerometer,  G  is  the  local 
value  of  gravity,  and  0  is  90°  minus  the  angle  the  sensitive  axis  makes  with  the 
local  vertical.  As  the  angle  with  the  vertical  becomes  larger  and  approaches  90°, 

6  becomes  0°  and  A  becomes  0.  A  90°  dip  is  all  that  is  required  to  calibrate  from 
zero  to  one  G.  For  this  experiment  one  need  only  calibrate  from  0.0001  to  0. 1  G, 
which  represents  an  angle  of  about  6°. 

The  small  angle  is  determined  very  accurately  with  the  help  of  a  double  sine  bar 

_  c 

in  an  area  where  the  noise  level  is  below  10  G.  The  linearity  and  repeatability  of 
the  accelerometer  were  measured  by  attaching  the  accelerometer  to  the  face  plate 
of  a  Leitz  dividing  head  and  computing  the  values  of  acceleration  from  the  sine  of 
the  angle.  (See  Figures  5  and  6.) 


3.  PREDICTION  OF  DRAG  COEFFICIENT 

The  sphere,  in  the  altitude  range  from  160  km  to  90  km,  falls  through  all  the 
recognized  regimes  of  rarefied  gas  flow,  being  in  free -molecular  flow  at  the  top 
and  in  laminar  boundary-layer  flow  at  the  bottom  of  this  range.  In  general  CD  will 
depend  upon  the  Mach  number  (or  molecular  speed  ratio),  Reynolds  number,  ratio 
of  sphere  temperature  to  ambient  temperature,  atmospheric  composition,  and  the 
'reflection  coefficients'  that  characterize  the  degree  to  which  molecules  leave  the 
sphere  surface  in  thermodynamic  equilibrium  with  it.  Fortunately,  not  all  these 
factors  are  important  in  every  range  of  altitude. 

To  make  specific  estimates  of  CD  vs_  altitude,  one  must  first  calculate  M,  Re, 
Kn,  and  Tw/Tq  vs  H,  using  best  prior  estimates  of  atmospheric  structure  at  the  time 
of  flight.  For  this  purpose,  the  ARDC  Model  Atmosphere,  1959  was  used.  Vis¬ 
cosities  and  mean  free-paths  were  estimated  from  a  simple  Sutherland  formula 
(see  Table  1,  page  11)  that  ignores  changes  in  atmospheric  composition. 


Figure  5.  Calibrating  Accelerometers  with  Double  Sine  Bar 


Figure  6.  Obtaining  Linearity  of  the  Accelerometers  with  a 
Laltx  Dividing  Head 
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According  to  these  estimates  the  sphere  was  in  free-molecular  flow  at  140  km  and 
above.  At  altitudes  where  drag  data  could  be  sensed  with  reasonable  accuracy, 

135  to  95  km,  the  sphere  was  in  the  transition  regime  and  the  boundary-layer  regime. 

3. 1  Near-Continuum  Flow 

For  altitudes  between  110  and  90  km  the  sphere  encountered  conditions  of  Mach 
number,  Reynolds  number,  and  T  /T  very  nearly  identical  to  those  investigated  by 

rj  W  O 

Aroesty  in  the  University  of  California  low-density  wind  tunnel.  His  results  are 

8  9 

substantially  corroborated  by  the  experiments  of  Ashkenas  and  Sreekanth.  While 
the  absolute  temperature  level  of  the  sphere  flight  was  higher  than  that  of  the  wind- 
tunnel  tests,  no  appreciable  deviations  from  thermally  perfect  gas  behavior  should 
have  appeared  in  either  case.  The  dependence  of  viscosity  upon  temperature,  which 
influences  the  conversion  between  Re^  and  Re2<  differs  appreciably  between  flight 
and  wind-tunnel  conditions;  and  it  has  been  assumed  that  Re„,  rather  than  Re  ,  is 
the  important  Reynolds  number  for  this  problem.  In  Figure  7  the  predicted  curves 
of  Cq  vs  Re2  are  essentially  reproductions  of  Aroesty's  data  for  Re2  >  22  at  M  ■  6. 
Re2  >  7  at  M  *  4,  and  Re2  >  6  at  M  =  2.  The  curve  for  M  »  3  is  interpolated. 

3. 2  Free-Molecular  Flow 

Above  140  km  the  Knudsen  number  exceeded  five,  and  free -molecular  flow  was 
assumed.  Although  no  drag  data  were  obtained  at  these  altitudes,  a  few  observations 
will  be  made  about  the  prediction  of  Cp  for  free -molecular  conditions,  since  they 
may  be  useful  in  planning  future  experiments  of  this  type  and  they  are  used  here  in 
an  interpolation  formula  for  transition  flow. 

The  basic  free -molecular  drag  theory  for  a  pure  gas  is  given  by  Schaaf  and 
Chambre'. 10  Upon  the  assumption  that  each  component  of  a  mixed  gas  acts  independ¬ 
ently ,  the  drag  coefficient  for  a  sphere  is  given  by 

cd=  ;nimicD(»i'»wi)/;nimi' 

where  n.  and  m.  are  the  number  density  and  particle  mass  of  the  ith  species. 

11  /  10 
each  species  the  authors  use  Equation  (8-6)  of  Schaaf  and  Chambre. 


(2) 

For 


DRAG  COEFFICIENT  C 


Figure  7.  Predicted  Curve#  of  vs  Re 
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where 


V  la  the  sphere  velocity,  relative  to  the  atmosphere;  Tw  is  the  sphere  surface 
temperature,  assumed  to  be  uniform;  and  T  is  the  kinetic  temperature  of  the 
atmosphere,  which  is  assumed  to  have  a  local  Maxwellian  distribution  of  molecular 
velocities.  T  is  assumed  to  be  the  same  for  all  molecular  species,  and  all  species 
are  supposed  to  be  electrically  neutral. 

The  principal  questionable  assumption  underlying  these  equations  is  that  each 
species  of  particle,  after  striking  the  sphere,  is  re -emitted  in  a  Maxwellian 
effusion  at  the  temperature  and  macroscopic  velocity  of  the  sphere,  but  with  un¬ 
changed  molecular  or  atomic  identity.  The  accuracy  of  this  assumption  is  expected 
to  be  quite  good  for  surface  materials  such  as  the  fabric  used  on  the  present  sphere, 
which  is  microscopically  rough  but  chemically  noncatalytic  for  oxygen  recombination. 

An  important  implication  of  Eq.  (3)  is  that  becomes  quite  insensitive  to 
atmospheric  temperature  and  composition  if  the  flight  speed  is  such  that  all  the 
speed  ratios  Sj  are  >  3.  On  the  other  hand,  the  nature  of  the  reflection  or  re- 
emission  of  molecules  at  the  sphere  surface  can  be  reasonably  well  predicted  from 
laboratory  experiments  if  the  average  kinetic  energy  of  the  atmospheric  particles 
relative  to  the  sphere  is  less  than  about  1  ev.  For  the  altitude  range  from  100  to 
200  km.  a  sphere  speed  of  about  3000  m/sec  would  be  ideal  in  both  respects.  As 
examples  of  the  errors  in  Cp  that  result  from  uncertainties  about  temperature  and 
composition,  three  values  of  Cp  for  h  *  160  km  have  been  contrasted.  The  value 
Cp  *  3.  81  given  in  Table  1  comes  from  using  the  ARDC  1959  Model  Atmosphere  and 
treating  the  atmospheric  mixture  as  a  pure  gas  having  a  mean  molecular  weight. 
Using  the  NASA  1962  (Harris  and  Priester* 1)  atmospheric  model,  and  again  treating 
the  mixture  as  an  equivalent  pure  gas,  gives  Cp  *  3.  33.  Finally,  use  of  the  NASA 
1962  model  and  treating  Ng.  Oj,  and  O  separately  will  give  Cp  *  3.  33.  If  these 
same  calculations  are  repeated,  with  the  sphere  speed  doubled,  the  corresponding 
results  are  2.  54,  2.  43,  and  2.46,  reducing  the  total  spread  of  values  from  14  per¬ 
cent  to  4.  5  percent. 


3, 3  Translation  Flow 

Between  110  and  140  km  the  sphere  was  in  transition  {low,  for  which  very  little 
directly  pertinent  theory  or  experimental  data  exist.  The  nearly-free -molecular 
flow  theories  do  not  seem  to  be  successful  as  they  stand,  *  and  the  authors  resort 
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*See  Maslach  and  Schaaf  for  a  critical  review. 
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here  to  a  simple  interpolation  formula  between  Aroeaty's  data  and  free-molecular 

theory.  The  choice  of  formula  ia  guided  indirectly  by  the  recent  cylinder  drag  data 
12 

of  Maalach  and  Schaaf,  the  very  high  Mach  number  aphere  drag  data  of  Masson, 
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Morria  and  Bloxsom,  and  general  intuitive  considerations  of  smoothness.  For 
Re 2  <  4  the  authors  take 

CD(s,  sw,  Re2)  -  CD  (s,sw)  [1  -  Re2(0.040  +  0.032  a"4)]  .  (4) 

F»M. 

It  is  felt  that  this  may  apply  fairly  well  over  a  range  of  values  of  Tw/T  and  with 
little  dependence  upon  atmospheric  composition,  except  for  the  influence  of  these 

quantities  uu  Cn 

UF.M. 

The  initial  estimates  of  Cjy  denoted  C^°\  and  based  on  the  ARDC  1959  Model 
Atmosphere,  are  tabulated  in  Table  1. 

4.  EXPERIMENTAL  RESULTS 
4. 1  Rocket  and  Sphere  Trajectory 

The  Aerobee  150  for  this  experiment  was  launched  from  the  Eglin  Gulf  Test 
Range  on  Santa  Rosa  Island,  Florida,  on  7  December  1961,  at  23:15:00.898  GMT 
(Fastax).  Sphere  ejection  occurred  80. 7  sec  later  at  285,  000  ft  and  the  sphere  then 
coasted  to  a  peak  at  250. 7  sec  at  732, 000  ft.  Rocket  vehicle  peak  was  at  254. 1  sec 
at  740,000  ft.  The  best  trajectory  for  the  rocket  was  obtained  by  the  MPS-19  radar 
beacon.  MATTS  data  provided  accurate  sphere  time -space  positioning  which  was  used 
in  the  density  calculations.  This  data  checks  with  optical  tracking  as  shown  in 
Figure  8. 

Sphere  velocity  determinations,  shown  in  Figure  9,  were  done  on  a  best-fit  basis 
by  using  local  values  of  acceleration  due  to  gravity  and  horizontal  velocity  determined 
from  values  at  the  peak  of  the  trajectory.  The  known  horizontal  velocity  was  290 
m/sec  and  was  assumed  constant  during  the  period  of  interest.  Any  variation  in 
horizontal  velocity,  resulting  from  some  external  force,  would  be  registered  by  the 
accelerometers  mounted  within  the  sphere.  The  maximum  measured  acceleration 
during  the  period  of  interest  was  0. 1  G  for  about  1  sec  or  a  velocity  variation  of 
0.98  m/sec,  which  is  an  insignificant  contribution  to  the  previously  determined 
velocities. 


(Xi)  aanxinv 


Figure  8.  Comparison  of  Sphere  Altitude  vs  Time  (MATTS  and  Optical  Data) 
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Figure  9.  Sphere  Velocity  va  Time 

4. 2  Accelerometer  Measurement 

Figure  10  ia  the  telemetry  record  tor  altitude  130  to  95  km,  representing 
approximately  24  sec  of  the  descent  portion  of  this  flight.  These  are  the  measured 
drag  forces  used  in  calculating  the  density  profile  for  the  experiment.  The  descent 
portion  of  the  flight  is  selected  because  the  sphere  platform  is  more  stable.  The 
ascent  portion  of  interest  occurs  immediately  after  ejection  and  time  is  required  for 
structural  rigidity  and  transitional  accelerations  to  damp  out.  A  wider  range  of 
more  accurate  measurements  was  possible  on  descent  for  this  flight. 

The  telemetry  record  represents  the  measured  accelerations  along  the  X,  Y,  and 
Z  axes  of  the  sphere,  shown  as  A?,  A1Q,  and  A^  reapectively.  These  axes  are  independ¬ 
ent  of  the  fixed  earth  axes  used  to  determine  altitude  and  velocity.  The  accelerometer 
measurements  were  telemetered  in  a  nonlinear  fashion  with  the  greatest  sensitivity 
at  the  low  values  of  acceleration.  This  increases  the  dynamic  range  of  the  instrument, 
increasing  the  maximum  altitude  for  the  experiment.  As  can  be  seen  from  this 
telemetry  record,  the  altitude  span  from  130  to  95  km  is  the  extent  of  sensitivity. 

This  is  not  an  inherent  limit  but  rather  represents  a  particular  design.  At  95  km  the 
readings  are  saturated  and  at  130  km  the  telemetered  data  uncertainty  is  large.  The 
commutated  telemetry  aignol  at  the  top  of  Figure  10  includes,  among  other  informa¬ 
tion,  signals  from  three  additional  accelerometers.  These  accelerometers  were 
included  in  this  experiment  as  back-up  Instruments  and  were  also  aligned  to  measure 
accelerations  along  the  X,  Y,  and  Z  axes.  A  less  sensitive  input  voltage  subcarrier 
channel  was  used  for  the  commutated  signal.  All  six  units  performed  correctly  and 
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Figure  11.  Uncorrected  Values  taken  directly 
from  Figure  10 

a  useful  check  was  provided.  This  discounts  any  possibility  of  erroneous  readings  due 
to  instrument  malfunctions  and  certifies  the  curves  of  Figure  11  without  explaining 
any  discrepancies  between  measured  and  predicted  values. 

Figures  11  and  12  are  plots  of  measured  accelerations  in  G  units  vs  altitude. 
Curves  A?,  A1Q,  and  AJ4  represent  the  X,  Y,  and  Z  axes  respectively  of  the  spinning 
sphere.  Figure  11  presents  the  uncorrected  values  taken  from  the  telemetry  record. 
The  only  difference  is  that  of  plotting  accelerations  linearly  instead  of  in  a  nonlinear 
function  as  displayed  by  the  telemetry  record.  Figure  12  shows  the  corrected 
values  of  drag  acceleration  obtained  from  Figure  11.  The  question  of  major  im¬ 
portance  is  *What  accelerations  need  to  be  subtracted  from  the  curves  of  Figure  11 
to  obtain  true  aerodynamic  drag?  "  Figure  13  shows  the  vector  sum  or  total  drag 
acceleration  from  Figure  12. 


Figure  12. 
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Figure  13.  Sum  Total  Measured  Drag 
Acceleration  in  G  Units  vs  Altitude  in  Kilo* 
meters 


4. 3  System  of  Forces 

For  a  first  approximation  the  only  forces  experienced  within  a  rigid  spinning 
sphere  in  free  fall  are  drag  and  centrifugal  forces.  The  drag  is  related  to  velocity 
with  respect  to  earth  coordinates,  and  centrifugal  acceleration  arises  from  the  spin 
of  the  sphere.  It  is  not  necessary  to  relate  these  independent  sets  of  axes  if  the 
spin  axis  of  the  sphere  remains  arbitrarily  constant  with  respect  to  its  X,  Y,  Z 
system.  In  this  case  the  centrifugal  acceleration  is  determined  by  the  normal  dis¬ 
tance  to  the  spin  axis  R  and  the  spin  rate  w.  If  the  accelerometer  coincides  with 

2 

the  spin  axis,  then  this  spurious  acceleration  will  be  zero.  In  practice  ac*u>  R  is 
minimized.  The  drag  occurs  opposite  the  velocity  vector  of  the  sphere  with  respect 
to  the  ambient  air.  For  high  velocities  in  still  air  this  is  the  direction  of  fall. 

Because  the  sphere  is  spinning  about  an  arbitrary  axis,  the  drag  generally  appears 
on  all  accelerometers  as  a  sine  function.  When  the  sphere  is  entering  the  atmosphere 
the  magnitude  is  always  increasing.  It  is  relatively  easy  to  distinguish  between  the 
quantities  a  ■  ajj  +  a£.  This  enables  one  to  subtract  the  error  ac  from  Figure  11 
and  obtain  true  drag  in  Figure  12. 

The  system  of  forces  acting  is  considerably  more  complicated  than  assumed 
above,  resulting  in  accelerations  that  cannot  always  be  disregarded.  The  sphere 
may  be  considered  a  solid,  having  a  rotation.  The  effects  that  play  an  important 
part  in  this  regard  are:  (1)  the  Poisson  effect,  (2)  the  Magnus  effect,  and  (3)  the 
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gyroscopic  effect.  Consider  the  case  of  a  spin  axis  having  a  component  normal  to  the 
direction  of  flow.  For  a  right-hand  rotation  the  Poisson  effect  results  in  a  devia¬ 
tion  to  the  right,  the  Magnus  effect  results  in  a  deviation  to  the  left,  and  the 
gyroscopic  effect  may  result  in  a  deviation  in  either  direction.  For  small  spin 
rates  the  gyroscopic  effect  will  predominate.  To  calculate  the  magnitudes  of  these 
unwanted  inputs  to  the  accelerometer,  one  must  know  the  orientation  of  the  spin 
axis.  It  is  possible  to  show  that  the  sum  effect  for  most  of  the  downleg  period  of 
interest  is  small. 

4. 4  Drag  Acceleration  Results 

To  consider  the  practical  circumstances,  the  sphere  did  not  act  as  a  rigid  body. 
While  the  sphere  instrument  package  was  nearly  20  pounds,  about  two-thirds  of 
this  weight  was  concentrated  in  the  instrument  package  located  at  the  center  and 
supported  by  two  inflatable  6-inch  struts  fastened  to  the  sphere  at  opposite  points 
on  the  equator.  This  method  of  suspension  is  easier  to  package  and,  after  ejection, 
should  inflate  easily  with  a  minimum  of  disturbance  to  the  structure.  It  appears  that 
on  a  zero  G  trajectory  (above  140  km)  the  system  damps  out  sufficiently  on  only  one 
axis— that  which  is  most  rigidly  supported— along  the  strut.  This  is  the  Z  axis  of  A^. 

Accelerometer  A14  has  a  constant  error  corresponding  to  the  centrifugal  accele¬ 
ration,  a  result  of  displacement  from  the  spin  axis.  Neglecting  other  forces,  a 
constant  was  subtracted  from  the  values  of  Figure  11  to  obtain  Figure  12  values  for 
this  accelerometer.  This  may  also  be  considered  as  subtracting  the  acceleration 
function  experienced  above  140  km  where  drag  is  known  to  be  zero  for  the  entire 
downleg  portion.  This  function  is  the  zero  reference  for  drag  purposes. 

The  zero  reference  or  error  to  be  subtracted  from  A10  is  a  periodic  function. 

This  apparently  arises  from  harmonic  motion  of  the  strut  in  the  cross-axis  plane, 
a  result  of  nonrigidity.  The  Y  axis  zero  reference  is  subtracted  from  the  curve  of 
Figure  11  to  obtain  Figure  12.  Since  the  aerodynamic  drag  has  a  frequency  of  about 
1/3  cps  and  the  zero  reference  has  a  much  longer  period,  this  can  be  done. 

Accelerometer  A?  experiences  the  least  amount  of  drag.  The  zero  reference  is 
not  sufficient  when  extrapolated  below  120  km,  and  a  new  zero  reference  for  95  to  120 
km  is  indicated  in  Figure  11.  This  represents  the  values  subtracted  to  obtain  Figure 
12.  The  basis  for  determining  the  zero  level  is  the  practical  consideration  that  the 
drag  which  remains  must  have  the  same  frequency  component,  1/3  cps.  This  amounts 
to  a  graphical  Fourier  analysis,  removing  the  undesirable  frequency  components  and 
leaving  the  aerodynamic  drag  frequency. 

Figure  13  represents  the  vector  sum  of  the  X,  Y,  and  Z  results  of  Figure  12,  drag 
in  G  vs  altitude  in  kilometers.  Since  CD  and  sphere  velocity  have  been  determined 
it  is  now  possible  to  obtain  the  density  profile  as  determined  by  Bq.(l).  Figure  14 
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shows  the  calculated  density  in  kilograms  per  cubic  meter  vs  altitude  in  kilometers, 
and  Figure  IS  compares  results  with  the  ARDC  Model  Atmospheres  of  1956,  1959, 
and  1961. 

5.  ERRORS 

The  sources  of  error  in  the  calculated  density  are:  drag  coefficient  data,  neglect 
of  winds,  calculation  of  the  trajectory,  and  in  measured  values  of  A,  m,  and  a^. 

Drag  acceleration  a^  is  the  maximum  source  of  error  at  present.  The  uncertainty 
in  Sq  arises  not  from  the  measuring  instrument  but  from  the  interpretation  of  tele¬ 
metered  data  to  obtain  true  drag  acceleration.  This  is  particularly  difficult  in  the 
case  of  a  nonrigid  sphere.  Experimental  uncertainties  are  given  in  Table  1. 

5.1  Drag  Coefficient  Data 

If  the  assumptions  made  in  the  earlier  discussion  regarding  the  reflection  pro¬ 
cess  in  free-molecule  flow  are  correct  above  130  km  (and  the  authors  feel  for  many 
reasons  that  they  are),  the  predicted  CQ  values  in  that  regime  should  have  no  more 
than  2  or  3  percent  probable  error  (due  possibly  to  uncertainties  in  surface  temper¬ 
ature).  The  estimated  probable  error  of  Cp  values,  where  they  actually  overlap 
Aroesty's  data  (below  110  km),  are  approximately  ±  2  percent.  There  are  no  firm 
guide  lines  for  estimating  error  in  the  interpolation  regime  between  free-molecule 
flow  and  Aroesty's  data.  Nonetheless  the  method  of  establishing  the  interpolation 
curves,  their  generally  reasonable  appearance,  and  the  fact  that  they  span  only  a 
relatively  narrow  gap  between  well-known  anchor  points  give  confidence  that  pre¬ 
dictions  made  from  these  curves  will  also  have  probable  errors  of  only  a  few  percent. 
It  would  take  a  major  distortion  of  any  of  these  curves  to  move  the  probable  error 
up  or  down  by  5  percent  of  Cq. 

5.2  Neglect  of  Wind 

Motion  of  ambient  air  with  respect  to  the  sphere  increases  the  effective  drag 
velocity  or  the  motion  of  the  sphere  with  respect  to  the  ambient  atmosphere.  Changes 
in  wind  speed  of  100  m/sec  have  been  observed  from  chemical  releases  within  a 
height  interval  of  10  km.  These  winds  are  mainly  horizontal  and  no  vertical  winds 
are  thought  to  be  significant  in  these  altitude  regions.  Neglecting  a  vertical  wind 
component  causes  a  much  larger  error  than  neglecting  a  horizontal  component  of  the 
same  magnitude.  However  the  nature  of  these  winds  is  such  that  a  variation  with 
altitude  usually  exists  and  wind  shears  are  often  encountered. 14  As  seen  in  Figure  13 
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the  Aerodynamic  drag  curve  ia  not  a  smooth  fit  of  any  exponential  curve.  Certain 
errors  are  introduced  in  the  ambient  velocity  especially  at  the  lower  altitudes  of 
the  flight  regime.  These  errors  are  included  in  the  total  error  estimates  of  Table  2. 


TABLE  2. 


Altitude 

(MATTS) 

km 

CD 

Velocity 

(MATTS) 

m/sec 

ad 

Total 

Drag 

ft 

Density 

kg/m3 

% 

Error 

Pred. 

CD 

% 

Error 

Neglect 

of 

Winds 

% 

Error 

in 

Traj. 

% 

Error  in 
Meas.  of 

aD 

Net  Prob. 
Error  in 
Calculated 
Density 

100 

1.38 

1591 

47 

5.94X10-7 

±2 

2 

>1 

±  6 

±  8% 

105 

1.66 

1561 

17.5 

1.88X10"7 

±2 

2 

>1 

±14 

±17% 

110 

2.05 

1530 

8.5 

8. 40x10"® 

±2 

2 

>1 

±  3 

±  4% 

115 

2.32 

1494 

4.5 

4.02X10-8 

5 

1 

>1 

±  4 

±  8% 

120 

2.49 

1454 

2.3 

2.09X10*8 

5 

1 

>1 

±  6 

±  8% 

125 

2.64 

1411 

1.2 

1.04X10"8 

5 

1 

>1 

±18 

±20% 

130 

2.78 

1366 

.6 

5.05X10'9 

3 

1 

>1 

±25 

±30% 

135 

2.90 

1326 

.3 

3.60X10-9 

3 

1 

>1 

±49 

±50% 

5. 3  Calculation  at  Trajectory 

Errors  in  the  calculation  of  the  trajectory  include  errors  in  altitude  that  affect 
the  position  of  the  data  on  Figure  13  and  errors  because  of  velocity  that  enter  into  the 
equation.  Position  uncertainty  of  the  sphere  during  descent  through  the  period  of 
interest  was  less  than  100  m.  Velocities  at  any  point  on  the  trajectory  have  two 
main  sources  of  error:  the  horizontal  velocity  at  peak  and  the  error  due  to  monitoring 
of  external  accelerations.  The  horizontal  component  was  determined  at  peak  to  be 
290  1  20  m/sec.  At  the  time  of  interest  the  vertical  velocity  is  large  and  the  con¬ 
tribution  of  the  horizontal  component  to  the  vector  sum  is  6  percent  or  less.  Thus 
the  total  error  in  velocity,  as  determined  by  the  earth  referenced  system,  is  less 
than  1  percent. 


5. 4  Drag  Acceleration 


Errors  in  A  and  m  can  be  made  negligibly  small  with  accurate  weighing  and 
good  dimensions  of  the  sphere.  Errors  associated  with  measurement  and  telemetry 
of  acceleration  are  a  function  of  magnitude  due  to  the  nonlinear  method  of  deduction. 
The  present  limit  on  maidmum  altitude  is  determined  by  the  telemetry  input  range. 
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At  the  appropriate  range  this  error  is  leas  than  1  percent.  As  greater  sensitivity 
is  desired,  then  the  methods  of  calibration  degrade  this  figure.  At  the  present  state 
of  the  art  these  results  could  be  extended  in  steps  to  altitudes  of  160  km.  Con¬ 
version  of  telemetered  acceleration  to  drag  values  accounts  for  the  unnecessarily 
large  error  (17  percent)  at  lower  altitudes.  This  is  due  to  the  nonrigidity  of  the 
sphere.  Consequently  the  spurious  accelerations  are  larger  and  an  uncertainty 
exists  as  to  their  exact  function  vs  time.  Determination  of  absolute  error  introduced 
in  aD  requires  knowledge  of  the  rigidity  of  the  sphere.  This  is  most  critical  at 
higher  values  of  a^.  At  115  km,  6  percent  over-all  uncertainty  in  the  experiment 
gives  a  picture  of  the  confidence  in  measurement.  Total  error  estimates  are  given 
in  Table  2. 

6.  ITERATIVE  IMPROVEMENTS  OF  Cd  ESTIMATES  AFTER  THE  FLIGHT 

From  estimates  of  vs  altitude,  and  with  the  measured  velocity  and  drag 

force  data,  a  first  estimate  of  the  density  at  the  specific  time  and  location  of  the 
sphere  flight  is  derived.  In  the  present  case,  the  densities  so  derived — referred  to 
as  p^ — were  substantially  larger  than  those  suggested  in  the  ARDC  Model  Atmosphere, 
1959.  This  fact,  together  with  the  evidence  of  the  importance  of  solar-cycle  and 
diurnal  variations  in  atmospheric  properties  above  100  kmf  makes  it  necessary  to  use 
the  sphere  drag  data  to  the  fullest  extent  possible  to  obtain  self-consistency  between 
the  properties  that  were  used  to  determine  and  those  deduced  from  it  and  the 
flight  data. 

The  obvious  iteration  procedure  is  to  use  p^(h)  and  an  assumed  value  of  RT  at 
the  highest  altitude  for  which  density  data  were  obtained  (135  km)  and  then  find 
(RT)^  from  the  combined  equations  of  hydrostatic  equilibrium  and  of  state.  Thus 

rh 

p(1)(RT)(1)  -  (RT)135  -  J  pgdt.  (5) 

lwd 

This  procedure  is  arbitrary  in  its  choice  of  RT  at  h  *  135  km  and  in  the  decision  to 
use  only  an  arbitrary  RT,  multiplying  it  by  a  known  Pjjg ,  rather  than  assuming  an 
arbitrary  p^-jj.  Equation  (5)  has  been  integrated  numerically,  approximating  the 
graph  of  p(D  g  vs  h  by  straigh-line  segments  for  each  Ah  =  5  km,  on  a  plot  of  log 
[p(1)g]  vs  h.  This  leads  to  the  formula 


*See,  for  example,  Harris  and  Priester.  ** 
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Several  calculations  were  made  to  assess  the  sensitivity  of  derived  estimates 

of  RT  to  the  assumed  RT  value  at  135  km  and  to  the  scatter  of  the  measured  p^ 

5  2  2 

values  around  a  smooth  curve.  In  particular  (RT)135  »  2. 23  X  10  km  /sec  was 
tried  from  the  ARDC  1959  Model  and  (RT)135  ■  1.48  X  105  km2/sec2  was  tried 
from  the  NASA  1962  Model  of  Harris  and  Priester.  * 1  The  latter  value  corresponds 
to  Harris  and  Priester's  solar  activity  parameter  s  «  100  (late  1961)  and  to  a  local 
flight  time  of  19:00  hours.  Figure  16  shows  the  results  of  these  integrations.  As 
expected  the  curves  converge  and  become  essentially  independent  of  (RT)133  (and 
of  for  the  uppermost  two  or  three  points,  which  were  near  the  limit  of  experi¬ 
mental  sensitivity)  where  p^  (h)  has  become  much  greater  than  Thus  the 

two  curves  start  at  (RT>135  (ARDC)  ■  1.  51  (RT)13g  (NASA),  but  are  only  7  percent 
apart  at  h  »  115  km  and  less  than  1  percent  apart  at  h  *  100  km. 

For  the  calculation  of  C^  ^  the  authors  chose  to  use  the  NASA  (RT)133>  assuming 
this  more  recent  prediction  to  be  more  applicable  than  the  ARDC  1959  model.  This 
resulted  in  a  sizable  increase  in  the  Mach  numbers  and  nearly  doubled  Re2  at  100  km 
while  slightly  reducing  it  at  135.  The  net  effect  on  CD  was  a  fairly  uniform  decrease, 
ranging  from  4.  8  to  8.  9  percent.  These  values  were  then  used  to  obtain  second 

estimates  of  density,  p‘2\  which  were  5  to  9  percent  higher  than  The  hydro- 

fit  (2) 

static  equation  was  again  integrated,  using  p'  ,  with  the  result  that  (RT)'  7  dif¬ 
fered  at  most  by  2  percent  from  (RT)^\ 

Finally  a  second  iteration  was  performed  to  test  convergence,  and  the  resulting 
(3)  to) 

density  estimates,  p'  differed  by  less  than  1  percent  from  p'  \  The  final 

values  cited  as  the  experimental  results  of  this  flight  are  p^  and  RT^2\  based  on 

graphically  smoothed  values  of  p^\  Figure  17  shows  this  smoothed  curve,  together 

with  the  original  point  values  and  their  probable  errors. 
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Figure  16.  First  Estimate  b  of  RT  based  on  Figure  17.  Final  Values  of  Experimental  Results 

Sphere  Drag  Density  Data  p'U  and  Assumed  RT 
at  h  *  135  kmlpU)  data  smoothed] 
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